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Abstract 

In this work, we present an alternate set of basis functions, each defined 
over a pair of planar triangular patches, for the method of moments solu- 
tion of electromagnetic scattering and radiation problems associated with 
arbitrarily-shaped, closed, conducting surfaces. The present basis functions 
are point-wise orthogonal to the pulse basis functions previously defined. 
The prime motivation to develop the present set of basis functions is to uti- 
lize them for the electromagnetic solution of dielectric bodies using a surface 
integral equation formulation which involves both electric and magnetic cur- 
rents. However, in the present work, only the conducting body solution is 
presented and compared with other data. 
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1 Introduction 


In a recent paper, we have introduced a new set of pulse-type basis functions, 
each defined over a pair of planar triangular patches, to calculate the electro- 
magnetic scattering/radiation associated with three-dimensional, arbitrarily- 
shaped conducting bodies [1]. The primary motivation for this work is to 
develop an efficient and well-conditioned method of moments (MoM) [2] so- 
lution for dielectric material bodies via a surface integral equation (SIE) ap- 
proach [3]. It may be noted that since the SIE approach involves both electric 
and magnetic currents as unknowns in the MoM formulation, it is necessary 
to define two mutually orthogonal sets of basis functions to generate a well- 
conditioned moment matrix. In our previous paper, we have presented one 
set of basis functions and their suitability to calculate electromagnetic scat- 
tering from arbitrarily-shaped conducting bodies. In this work, we present 
the alternate set of basis functions and the numerical solution scheme to 
solve the same class of problems. We emphasize here that the mathematical 
equations appearing in the dielectric body SIE formulation are very similar 
to the equations we encounter in the PEC body problem, and, therefore, 
solving the PEC case provides confidence in the solution methodology. 


2 Description of the Problem 

Let S denote the surface of an arbitrarily-shaped, perfectly conducting, closed 
body illuminated by an incident electromagnetic plane wave E' as shown in 
Figure 1. Using the equivalence principle, potential theory, and the free- 
space Green’s function [4], the electric field integral equation (EFIE) is given 
by 

\juA + V$]tan = E l tan ( 1 ) 

where the subscript “tan” refers to the tangential component. In (1), 
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e and /i are the permittivity and permeability of the surrounding medium, 
k is the wave number and r and r' represent the position vectors to ob- 
servation and source points, respectively, from a global coordinate origin. 
The unknown surface current Js is related to the charge density qs by the 
continuity equation, given by 

V • J s = ~ju qs • (6) 

For the numerical solution of (1), we apply the method of moments for- 
mulation using planar triangular patch modeling and the basis functions as 
described in the following: 


3 Description of Basis Functions 


Let T+ and T~ represent two triangles connected to the edge n of the triangu- 
lated surface model as shown in Figure 2. We define two mutually orthogonal 
vector basis functions associated with the n th edge as 


fn ( r ) 
and 
9 n (r) 


n ± x £, re S n , 
0, otherwise 


£, re S nt 

0, otherwise 


( 7 ) 
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where S n represents the region obtained by connecting the mid-points of the 
free edges to the centroids of triangles T 7 f and to the nodes of edge n. This 
area is shown shaded in Figure 2. Also, £ and h ± represent the unit vector 
along the n th edge and the unit vector normal to the plane of the triangle 
T±, respectively. Note that the basis functions defined in (7,8) are actually 
pulse functions, point-wise orthogonal to each other, defined over the region 
S n . The electromagnetic solution procedure using the basis functions f n is 
described in [1], and in this work we present a similar solution using the basis 
functions g n . 

It may also be noted here that the basis functions g n are less versatile than 
the basis functions f n . Unlike the basis functions f n , which are applicable 
to both open and closed bodies as demonstrated in [1] , the basis functions g n 
are applicable to closed bodies only. The main reason for this restriction is 
that the functions g n are defined parallel to the edges of the planar triangular 
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patches and represent the tangential component of the current. It is well- 
known that the tangential component of the surface current is undefined for 
an open surface at the boundary. However, our primary motivation to use 
these functions is their applicability in the solution of the dielectric body 
problem, which is always posed for a closed body. Hence, the basis functions 
g n do not hinder our purpose. 


4 Numerical Solution Procedure 

As a first step, we consider the testing procedure. Consider the m th interior 
edge, associated with triangles T±. We integrate the vector component of 
(1) along the m th edge to obtain 

[ juA(r) • dt+ I V<h(r)«d£ = f E\r)*dt (9) 

which may be re-written as 

L jcuA(r).d £+ = f E\r).di (10) 

for n = 1,2,3,- • -,N, where n represents the source charge region and N 
represents the total number of edges in the triangulation scheme. Note that 
in (10), <h™ 1 and <h™ 2 represent the scalar potentials evaluated at the nodes 
connected to the m th edge. Further, the integrals appearing in (10) may be 
easily evaluated using any accurate numerical integration algorithm such as 
the 1-point, 2-point or 4-point trapezoidal rule. In this work, we choose to 
use a 2-point method. 

Next, we consider the expansion procedure. Using the basis functions g n 
defined in (8), we approximate the unknown current J as 

N 

J = J2 l n 9n ■ (H) 

n = 1 

This is followed by a substitution of the current expansion (11) into (10), 
yielding a, N x N system of linear equations which may be written in matrix 
form as ZI = V, where Z = [Z mn ] is an JV x JV matrix and I = [ I n ] and 
V = [V m ] are column vectors of length N. The elements of the Z and V are 
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given by 
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and r m i and r m2 are the position vectors to the nodes ml and m2, respec- 
tively, connected to the m th edge. 

The numerical evaluation of the vector potential, shown in (14), is straight- 
forward and may be accomplished by the procedure described in [5]. How- 
ever, the numerical evaluation of the scalar potential terms described in (15) 
and (16) may be carried out as follows: 

Let us define the unknown charge density q$ in (3) as 

N n 

qs = a i p i ( 20 ) 

i= 1 

where N n represents the total number of nodes (vertices) in the model and 


p i(r) 


1 , re Si 

0, otherwise. 


( 21 ) 


In (21), Si is the i th charge patch, formed by connecting the centers of the 
edges and the centroids of the triangles associated with the node i, as shown 
by the shaded area in Figure 3. 
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Now, let us consider the i th charge patch. Using (6), the well-known 
Divergence theorem and simple vector calculus, we have 



f 

JSi -ju 


LO 


f J • fig <M 

JCi 
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where C % is the contour bounding the charge patch S), and hi is the unit 
vector normal to the contour Ci in the plane of the triangle containing the 
contour segment. By considering the left hand side of (22), we have 




cp Asi 


(23) 


where Asi represents the area of the charge patch S'*. Considering the right 
hand side of (22), we have 
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where £ tJ , r^, and hf 3 represent the unit vector along the j th edge connected 
to node i, the vector from the center of the j th edge to the centroid of 
the ± triangle connected to the edge, and the unit vector normal to the 
plane of the ± triangle associated with the j th edge connected to the node 
i, respectively. Also, Ek represents the total number of edges connected to 
node i. Using (23) and (24), we have 

j e k 

= — r- J2 Aj[4- • ( r p + x A tj + ( r ii x Kj)\ ■ ( 25 ) 

-Si j = l 


Thus, we can write the scalar potential terms in (15) and (16) as 
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where r m i , r m2 , r ni , and r n2 are the position vectors to nodes ml and m2 
connected to the edge rn, and nodes nl and n2 connected to the edge n, 
respectively. Further, S n i and S n2 are the charge patches associated with the 
nodes nl and n2, respectively, as depicted in Figure 3. The integrals in (26) 
and (27) may be evaluated with the procedures described in [5]. 

Finally, once the matrices Z and V are determined, we may easily solve 
the system of linear equations to obtain I. 


5 Numerical Results 

In this section, we present numerical results for a PEC sphere (diameter 
= 0.15A) and a cube (length = 0.15A), and we compare the results with 
the solution obtained using the procedure presented in [1], Also, for the 
case of the sphere, the results are compared with the exact solution. The 
sphere and the cube are modeled with 500 and 960 triangles, respectively. 
For both examples, the body is placed symmetrically at the center of the 
coordinate system and illuminated by an x-polarized plane wave traveling 
along the z-axis. The bistatic radar cross section is presented in Figures 4 
and 5. We note that the new basis function results compare well with the 
other methods. 


6 Conclusions 

In this work, we present a new set of basis functions for the method of 
moments solution of electromagnetic scattering by conducting bodies of ar- 
bitrary shape. The new basis functions are pulse basis functions defined over 
a pair of triangular patches and tangential to the common edge. The present 
set of basis functions along with the functions defined in our previous work 
[1] should prove helpful in obtaining a stable and well-conditioned solution 
to the material body problem. At present, the work is in progress to apply 
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these two sets of basis functions to material bodies, and the results will be 
reported in the near future. 
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Figure 1: Arbitrarily-shaped conducting body excited by an incident elec- 
tromagnetic plane wave. 
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Figure 2: Basis function description: For the n th edge, basis functions f and 
g are defined in the shaded area. 
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Figure 3: Charge patch around i th node. 
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Figure 4: Bistatic RCS of a sphere of diameter 0.15 A excited by a plane 
wave travelling in the -z direction. 
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Figure 5: Bistatic RCS of a cube of length 0.15 A excited by a plane wave 
travelling in the -z direction. 
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